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1 Introduction 



Quantum field theory lies at the root of modern physics. After the success of the 
standard model in describing particle physics, one of the most pressing open question 
is how to derive an extended version of field theory which encompasses the quanti- 
zation of gravity. There are several attempts for this, among which string theory, 
loop gravity and noncommutative geometry are the best known. In each of these 
attempts one of the key problem is to relax the constraints that formulate quantum 
field theory on a particular space-time geometry. 

What is certainly more fundamental than geometry is topology and in particular 
discrete structures on finite sets such as the species of combinatorists [IJ. The most 
prominent such species in field theory is the species of Feynman graphs. They were 
introduced by Feynman to label quantum field perturbation theory and to automatize 
the computation of connected functions. Feynman graphs also became an essential 
tool in renormalization, the structure at the heart of quantum field theory. 

There are two general canonical operations on graphs namely the deletion or con- 
traction of edges. Accordingly perhaps the most important quantity to characterize 
a graph is its Tutte polynomial [2j, l3j . This polynomial obeys a simple recursion rule 
under these two basic operations. It exists in many different variations, for instance 
multivariate versions, with possible decorations at vertices. These polynomials have 
many apphcations, in particular to statistical physics. For recent reviews see IHElin]- 

In recent years the Tutte polynomial has been generalized to the category of 
ribbon graphs, where it goes under the name of the BoUobas-Riordan polynomial 
[ZllHllG]. Around the same time physicists have increasingly turned their attention to 
quantum field theory formulated on noncommutative spaces, in particular flat vector 
spaces equipped with the Moyal-Weyl product [H] . This type of quantum field theory 
is hereafter called NCQFT. It happens that perturbation theory for such NCQFT's 
is no longer labeled by ordinary graphs but by ribbon graphs, suggesting a possible 
connection to the work of Bollobas-Riordan. 

Quantum field perturbation theory can be expressed in several representations. 
The momentum representation is the most common in the text books. The direct 
space representation is closer to physical intuition. However it is the parametric rep- 
resentation which is the most elegant and compact one. In this representation, after 
the integration of internal position and/or momentum variables has been performed 
explicitly, the result is expressed in terms of the Symanzik polynomials. There is 
an extensive literature on these polynomials (see e.g. [TOl [11] for classical reviews). 
These polynomials only depend on the Schwinger parameters. Space time no longer 
enters explicitly into that representation except through its dimension which appears 
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simply parameter. 

This observation is crucial for several key applications in QFT which rely on 
dimensional interpolation. Dimensional regularization and renormalization was a 
crucial tool in the proof by 't Hooft and Veltmann that non-Abelian gauge theories 
are renormalizable [1^. The Wilson-Fisher e expansion |13] is our best theoretical 
tool to understand three dimensional phase transitions. Dimensional regularization 
is also used extensively in the works of Kreimer and Connes [HI [T5] which recast the 
recursive BPHZ forest formula of perturbative renormalization into a Hopf algebra 
structure and relate it to a new class of Riemann-Hilbert problems [16]. 

Following these works, renormalizability has further attracted considerable in- 
terest in the recent years as a pure mathematical structure. The renormalization 
group ambiguity reminds mathematicians of the Galois group ambiguity for roots 
of algebraic equations [IT]. Hence the motivations to study quantum field theory 
and renormalization come no longer solely from physics but also at least partly from 
number theory. 

The fact that the parametric representation is relatively independent of the de- 
tails of space time makes it also particularly appealing as a prototype for the tools 
we need in order to quantize gravity. The point of view of loop gravity is essentially 
based on the diffeomorphism invariance of general relativity. In the spin foam or 
group field theory formalism amplitudes are expressed as discrete sums associated to 
combinatoric structures which generalize Feynman graphs. They are in fact gener- 
alizations of ribbon graphs. To extend the parametric representation and eventually 
the theory of renormalization to this context is a major challenge, in which some 
preliminary steps have been performed [18] . 

In this paper we uncover the relationship between universal polynomials of the 
Tutte and BoUobas-Riordan type and the parametric representation in quantum field 
theory. The Symanzik polynomials that appear in ordinary commutative QFT are 
particular multivariate versions of Tutte polynomials. The relation between BoUobas- 
Riordan polynomials and the non commutative analogs of the Symanzik polynomials 
uncovered in |19l|20l|2l] is new. This establishes a relation between NCQFT, combi- 
natorics and algebraic topology. Recently the relation between renormalization and 
topological polynomials was explored in [22|, and in [2S]- We intend also to explore 
in the future the relation between Feynman amplitudes and knot polynomials. 

The plan of this paper is as follows. In the next section we give a brief introduc- 
tion to graph theory and to Tutte-like polynomials. In the third section we derive the 
parametric representation of Feynman amplitudes of QFT and give a new method to 
compute the corresponding Symanzik polynomials. The deletion/contraction prop- 
erty (2.11) of these polynomials is certainly not entirely new [2H 125] • But our 
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method which starts from the phase-space representation of Feynman amphtudes 
is inspired by earher work on NCQFT [20l |2T] and introduces two main technical 
improvements. One is the use of Grassmann variables to exploit the quasi-Pfaffian 
structure of Feynman amplitudes. This quasi-Pfaffian structure was discovered in 
|20] in the context of NCQFT but to our knowledge was never applied to the simpler 
case of ordinary QFT. The second improvement is that we do not factor out as usual 
the delta functions expressing global momentum conservation, because this requires 
a noncanonical choice of a root for every connected graph. Instead we introduce 
an infrared regularization in the form of a small harmonic potential at each vertex 
which leads to more elegant and canonical formulas. The corresponding generalized 
Symanzik polynomials obey a transparent deletion/contraction relation which allows 
to identify them with particular multivariate Tutte polynomials. These polynomials 
are close but not identical to the polynomials of [1]; we show how they both derive 
from a more general " categorified" polynomial. The usual Symanzik polynomials are 
simply recovered as the leading terms when the small harmonic potentials tend to 
zero. 

For completeness we also include a more standard way to compute the Symanzik 
polynomials through x space representation and the tree matrix theorem. 

In the fourth section we introduce ribbon graphs and Bollobas-Riordan polynomi- 
als. In the fifth and last section we define the first and second Symanzik polynomials 
of NCQFT and relate them to the Bollobas-Riordan polynomials, using again the 
Pfaffian variables. Formulas for such polynomials were first sketched in (TH], but 
without proofs, nor relation to the Bollobas-Riordan polynomials. 

In a companion paper we shall discuss generalizations of the Tutte and Bollobas- 
Riordan polynomials that occur for non-translation invariant theories with propaga- 
tors based on the Mehler rather than the heat kernel. These theories appeared as the 
first examples of renormalizable NCQFT's fI7\ ES| [221 EI] and they are the most 
promising candidates for a fully non-perturbative construction of a field theory in 
four dimensions |32l ESj [Ml |35] . In this case the harmonic potentials on the vertices 
are no longer needed as the Mehler kernel already contains an harmonic potential 
for the propagators of the graphs. 

2 Tutte Polynomial 

2.1 Graph Theory, Notations 

A graph G is defined as a set of vertices V and of edges E together with an incidence 
relation between them. The number of vertices and edges in a graph will be noted 
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edge flag self-loop 

Figure 1: Basic building blocks of a graph 



also V and E for simplicity, since our context always prevents any confusion. Graph 
theorists and field theorists usually have different words for the same objects so a 
little dictionary may be in order. We shall mostly use in this review the graph 
theorists language. In subsection |2.4| we introduce also some enlarged notion of 
graphs, with decorations called flags which are attached to the vertices of the graph 
to treat the external variables of physicists, plus other decorations also attached to 
vertices called (harmonic) weights to regularize infrared divergences. Generalizations 
to ribbon graphs will be described in section |4} 

Edges in physics are called lines (or propagators). Edges which start and end at 
the same vertex are definitely allowed, and called (self)-loops in graph theory and 
tadpoles in physics. A proper graph, i.e. a graph G without such self-loops, together 
with an arrow orienting each edge, can be fully characterized through its incidence 
matrix e„e- It is the rectangular E hj V matrix with indices running over vertices 
and edges respectively, such that 

• eye is +1 if e starts at v, 

• e„e is -1 if e ends at v, 

• e„e is otherwise. 

It is also useful to introduce the absolute value rj^e = |e^e| These quantities can be 
then generalized to graphs with self-loops by defining ee„ = for any self-loop e and 
vertex v but rj^v = 2 for a self-loop attached at vertex v and rj^v = otherwise. The 
number of half-edges at a vertex v is called the degree of v in graph theory, noted 
d{v). Physicists usually call it the coordination number at v. A self-loop counts for 
2 in the degree of its vertex, so that d{v) = J2e Vev 

An edge whose removal increases (by one) the number of connected parts of the 
graph is called a bridge in graph theory and a one-particle-reducible line in physics. 

A forest is an acyclic graph and a tree is a connected forest. A cycle in graph 
theory is a connected subset of n edges and n vertices which cannot be disconnected 
by removing any edge. It is called a loop in field theory. 
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Physicists understood earlier than graph theorists that half-edges (also called 
flags in graph theory [36]) are more fundamental than edges. This is because they 
correspond to integrated fields through the rule of GauBian integration, which physi- 
cists call Wick's theorem. Feynman graphs form a category of graphs with external 
flags decorating the vertices. They occur with particular weights, in physics called 
amplitudes. These weights depend on the detail of the theory, for instance the space- 
time dimension. A quantum field theory can be viewed the generating functional for 
the species of such weighted Feynman graphs. In this paper we shall reserve the 
convenient word flag exclusively for the "external fields" decorations and always use 
the word half-edge for the " internal half-edges" . 

An edge which is neither a bridge nor a self-loop is called regular. We shall call 
semi-regular an edge which is not a self-loop, hence which joins two distinct vertices. 

There are two natural operations associated to an edge e of a graph G, pictured 
in Figure |2} 

• the deletion, which leads to a graph noted G — e, 

• the contraction, which leads to a graph noted G/e. If e is not a self-loop, 
it identifies the two vertices fiand V2 at the ends of e into a new vertex Vi2, 
attributing all the flags (half-edges) attached to Vi and V2 to V12, and then it 
removes e. If e is a self-loop, G/e is by definition the same as G — e. 

A subgraph G' of G is a subset of edges of G, together with the attached vertices. 
A spanning forest of G is an acyclic subgraph of G that contains all the vertices of 
G. If G is connected a spanning forest is in fact a tree of G and any such spanning 
tree has — 1 vertices. 

As explained in the introduction a topological graph polynomial is an algebraic 
or combinatoric object associated with a graph that is usually invariant under at 
least graph homeomorphism. It encodes information of the graph and so enables 
combinatoric and algebraic method to deal with graphs. 

The Tutte polynomial [2] is one of the most general polynomial to characterize 
a graph. It is defined under a simple rule through the deletion and contraction of 
edges. It can be generalized to the larger theory of matroids [57] . 

The original Tutte polynomial which is a function of two variables can be gener- 
alized in various ways to multi-variable polynomials which have many applications, 
in particular in statistical mechanics where it evaluates the Potts model on graphs 
[H ISl [6] . These applications shall not be reviewed here. 

We present first the two main equivalent definitions of the Tutte polynomial. 
One direct way is to specify its linear recursion form under contraction of regular 
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Figure 2: The contraction-deletion of a graph 
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edges (which are neither loops nor bridges), together with an evaluation on terminal 
forms solely made of bridges and self-loops. Another definition IS ClS db rank-nullity 
generating function. By induction these definitions can be proved equivalent. 



2.2 Tutte Polynomial 

The definition through a recursion relation is a reduction rule on edges together with 
an evaluation for the terminal forms. The Tutte polynomial may be defined by such a 
linear recursion relation under deleting and contracting regular edges. The terminal 
forms, i.e. the graphs without regular edges are forests (i.e. graphs made of bridges) 
decorated with an additional arbitrary number of self-loops at any vertex. The Tutte 
polynomial evaluated on these terminal forms simply counts separately the number 
of bridges and loops: 

Definition 2.1 (Deletion-Contraction) If G = {V,E) is a graph, and e is a reg- 
ular edge, then 

Tcix^y) =TG/eix,y) + TG-eix,y). (2.1) 
For a terminal form G with m bridges and n self-loops the polynomial is defined by 

TG{x,y)=x^y\ (2.2) 



It is not obvious that Definition 2.1 is a definition at all since the result might de- 
pend on the ordering in which different edges are suppressed through deletion/contraction, 
leading to a terminal form. The best proof that Tq is unique and well-defined is in 
fact through a second definition of the Tutte polynomial as a global sum over sub- 
graphs. It gives a concrete solution to the linear deletion/contraction recursion which 
is clearly independent on the order in which edges are suppressed: 

Definition 2.2 (Sum overs subsets) If G = {V,E) is a graph, then the Tutte 
polynomial of G, TG{x,y) has the following expansion: 

TG{x,y) = $^(x - lYi^)'^i^\y - 1)«(^), (2.3) 

AcE 

where r{A) = \V\ — k{A) is the rank of the subgraph A and n{A) = |y4| + k{A) — \ V\ 
is its nullity or cyclomatic number. In physicists language n{A) is the number of 
independent loops in A. 

Remark that r{A) is the number of edges in any spanning forest of A, and n{A) 
is the number of remaining edges in A when a spanning forest is suppressed, so it is 
the number of independent cycles in A. 
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Theorem 2.1 These two definitions are equivalent. 

One can show that the polynomial defined by the sum over subsets obeys the 
deletion-contraction recursion. One can also evaluate it directly and show that it 
coincides with the first definition on the terminal forms with only loops and bridges. 

There is a third definition of the Tutte polynomial through spanning trees (see 
eg [S]). This third definition involves ordering the edges of the graph. We think it 
may be also relevant in the context of field theory, in particular in relation with the 
ordered trees or forests formulas of constructive theory [3H1 EHl SO] , but this point of 
view will not be developed here. 

2.3 Multivariate Tutte polynomials 

Multivariate Tutte polynomials can also be defined through linear recursion or global 
formulas. 

The ordinary multivariate Tutte polynomial Zciq, {/?}) has a different variable 
/?e for each edge e, plus another variable q to count vertices. We also write it most of 
the time as ZG{q,l3) for simplicity. It is defined through a completely general linear 
deletion-contraction relation: 

Definition 2.3 (Deletion-Contraction) For any edge e (not necessarily regular) 



This relation together with the evaluation on terminal forms completely defines Zclq, (3), 
since the result is again independent of the order of suppression of edges. The termi- 
nal forms are graphs without edges, and with v vertices; for such graphs Zclq, j3) = q^ . 

We can also define ZG{q,P) as a sum over subsets of edges: 

Definition 2.4 (Sum over subsets) 



where k{A) is the number of connected components in the subgraph {V,A). 

One can prove as for the two variables Tutte polynomial that this definition is equiv- 
alent to the first. In [1] this multivariate polynomial is discussed in detail. 



Zciq, {P}) = PeZc/eiq, {P - {Pe}}) + ^G-e(g, {P ' {Pe}})- 



(2.4) 




(2.5) 
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To understand the relation between this multivariate and the ordinary Tutte 
polynomial with two variables we multiply Zq by we set /3e = ?/ — 1 and 

q= [x — — V) and get 

[g-^Za(g,/?)]|/3.=,-i,,=(-i)(.-i) = (x-1)'(^)-I^ITg(x,i/). (2.6) 
We consider also 



■'(''^^G(g,/5). (2.7) 



Taking the limit g — > that is retaining only the constant term in q we obtain a sum 
over maximally spanning subgraphs A, that is subgraphs with fc(A) = k{G\. 

Sg{(3)= Y1 (2.8) 

A maximally spanning E e^A 

If we now retain only the lowest degree of homogeneity in (3 we obtain a sum over 
maximally spanning graphs with lowest number of edges, ie maximally spanning 
acyclic graphs or spanning forests of G. 

Fg{P)= E Uf^e (2.9) 

maximally spanning forest of G e6JF 

Finally if we divide Fc{P) by YleeE /^e and change variables to = we obtain 
the "(Kirchoff-Tutte)-Symanzik" polynomial. This polynomial is usually defined for 
connected graphs, in which case the sum runs over spanning trees T of G. 

T spanning tree of G e^T 

This polynomial satisfies the deletion contraction-recursion 

UG{a) = Uc/eia) + aeUc-eia) (2.11) 
for any regular edge e, together with the terminal form evaluation 

UGia)= H a„ (2.12) 

e self— loop 



for any G solely made of self-loops and bridges. The deletion-contraction (2.11) can 
be extended to general edges if we define U for disconnected graphs as the product 
over the connected components of the corresponding f/'s and put the contraction of 
any self-loop to 0. 
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The polynomial U appears in a key computation of QFT, namely that of the 
parametric representation of the Feynman amplitude associated to the graph G. We 
give a proof of this fact based on a new Pfaffian representation of Feynman amplitudes 
together with harmonic weights at vertices so as to make the deletion/contraction 



rule (2.11 )-(2.12) particularly transparent. 

But to define the second (Kirchoff-Tutte)-Symanzik polynomial as well as to make 
the computation of the first Symanzik polynomial more canonical, we need first to 
enlarge slightly our category of graphs to include some decorations at the vertices. 

2.4 Decorated graphs 

Decorations are essential in physics to represent the concept of external variables, 
which are ultimately those connected to actual experiments and observations. 

Graphs with integers attached to each vertex and their corresponding multivariate 
polynomials WciaeyNy) have been considered in [13]. But to represent external 
variables we need to replace the integer Ny by a set of Ny disjoint object^ hereafter 



called flags (see subsection 2.1) 



Each flag is attached to a single vertex. A momentum variable pj in M*^ is associ- 
ated to each such flag. The incidence matrix can be extended to the flags, that is we 
define e/„ as +1 if the flag / is associated to the vertex v and otherwise. The total 
momentum incident to a subset S of the graph is then defined as Yliv&s^f'^'Pf- 
mark that this momentum is defined for subgraphs S which may contain connected 
components reduced to single vertices. For translation invariant QFT's, global mo- 
mentum conservation means that the condition = must be fulfilled. 

Similarly we attach to each vertex a number qv > ^ called the (harmonic) weight 
of the vertex. The total weight of a subgraph S is Ylv^^s I'"- 

The deletion/contraction relation is then extended to this category of graphs. 
The deletion is easy but the contraction is a bit non trivial. For a semi-regular edge 
joining vertices vi and V2 it collapses the two vertices into a single one vu, attaching 
to f 12 all half-edges of vi and V2- But it also attaches to vu the union of all the flags 
attached to Vi and V2, so that the total momentum incoming to f 12 is the sum of the 
momenta incoming to f 1 and to f 2. Finally the new weight of f 12 is the sum + 
of the weights of Vi and V2- 

These decorated graphs are the natural objects on which to define generalized 
Symanzik polynomials in field theory. 

Remaining for the moment in the context of graph theory we can define the 
second (Kirchoff-Tutte)-Symanzik polynomial for a connected graph as 



^In mathematics such a replacement is called a categorification of the integers Ny. 
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Definition 2.5 



Vcia^p) = -^^p^ -p^' JJ (2.13) 

v^v' 72 2— tree separating v and e^72 

where a two tree 7^ means a tree minus one edge, hence a forest with two disjoint 
connected components Gi and G2; the separation condition means that v and v' must 
belong one to Gi the other to G2- 

For any pair of distinct vertices v and v' we can build a canonical graph G{v, v') 
first by joining vertices v and v' in G with a new edge and then contracting that edge. 
This operation could be called the contraction of the pair of vertices v and v' . The 
following result goes back to Kirchhoff |41j. 

Proposition 2.1 The second Symanzik polynomial is a quadratic form in the total 
momenta p^ at each vertex, whose coefficients are the Ug(v,v') polynomials: 

VG(a,p) = -Pt,' Ug{v,v')- (2.14) 

Proof The graph G{v, v') has V — 1 vertices, hence its spanning trees have V — 2 
edges. They cannot make cycles in G because they would make cycles in G{v,v'). 
They are therefore two-trees in G, which must separate v and v', otherwise they 
would make a cycle in G{v,v'). □ 
On the submanifold of flag variables satisfying the momentum conservation con- 
dition pg = ^fPf = there is an alternate less symmetric definition of a similar 
polynomial: 

Definition 2.6 

7^2 2-trce e^T2 

where T2 is again a two-tree with two disjoint connected components Gi and 6*2 ■ 
Indeed this is an unambiguous definition. On the submanifold = we have 



PGi = — PG2! hence equation (2.15) does not depend of the choice of Gi rather than 
G2. 

Proposition 2.2 On the manifold of flag variables satisfying the momentum con- 
servation condition pg = J2fPf — ^''^^ ^^(a,^) = Vg^cx.p). 
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Proof We simply commute the sums over v,v' and 7^ in (2.13). For a given 7^ the 
condition that v and v' are separated allows to separate the with v E Gi from the 
p^i with v' G (j2; one gets therefore — \ 2pG^ ■PG2 which is nothing but Pq_^ or Pq^ on 
the manifold pc = 0. □ 



We shall give in subsection 3.4 a definition of generalized first and second Symanzik 
polynomials for any graph, connected or not from which Ug, Vq or Vq can be easily 
derived in certain limits. Before actually performing these computations we include 
a brief interlude on Grassmann representation of determinants and Pfaffians. The 
reader familiar with this topic can jump directly to the next section. 

2.5 Grassmann representations of determinants and Pfaffi- 
ans 

Independent Grassmann variables xii---iXn satisfy complete anticommutation rela- 
tions 

XiXj = -XjXi Vi,j (2.16) 

so that any function of these variables is a polynomial with highest degree one in 
each variable. The rules of Grassmann integrations are then simply 



dx = 0, Jxdx = l- (2.17) 

The determinant of any n by n matrix can be then expressed as a Grassmann 
GauBian integral over 2n independent Grassmann variables which it is convenient to 
name as ipi, . . . jipn, ipi, . . . yipn, although the bars have nothing yet at this stage to 
do with complex conjugation. The formula is 

det M = diJid^iC- ^^^^^^J . (2. 18) 

The Pfaffian Pf(A) of an antisymmetric matrix A is defined by 

detA=[Fi{A)f. (2.19) 
Proposition 2.3 We can express the Pfaffian as: 

Pf(A) = j dxi...dxne~^^<^^'^''^' = j dxi-.dxne'^^^-^^''^''^'. (2.20) 
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Proof Indeed we write 

detA= fjldiPidiJie-^^^'^^^^^'^^. (2.21) 

Performing the change of variables (which a posteriori justifies the complex notation) 

(2.22) 



whose Jacobian is z the new variables x and u are again independent Grassmann 
variables. Now a short computation using Aij — — -^ji gives 

detA = f Yldxidujie-^'<y>^"^^^^^'^'<^'^^'^'^'^^ 

i 

Yl dxiC ^^<^ Yl ^"^i^" ' (2.23) 

i i 

where we used that n = 2p has to be even and that a factor (— 1)*^ is generated when 
changing Yl - dxidui into J^- dxi Hi ^^i- Equation (2.23) shows why det A is a perfect 
square and proves (2.20). □ 

Lemma 2.1 The determinant of a matrix D + A where D is diagonal and A anti- 
symmetric has a "quasi-Pfaffian" representation 



det{D + A) = j Y[dXiduJie-^^^^^"'^'-^'<^^''^^^^^-^^^<i'^'^^^'^^. 



(2.24) 



Proof The proof consists in performing the change of variables (2.22) and canceling 
carefully the i factors. □ 



2.5.1 Tree-Matrix Theorem 

Let A be an n X n matrix such that 

n 

= Vj . (2.25) 

i=l 

Obviously det A = 0. The interesting quantities are eg the diagonal minors det A" 
obtained by deleting the i-th row and the i-th column in A. The "Kirchoff-Maxwell" 
matrix tree theorem expresses these minors as sums over trees: 
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Theorem 2.2 (Tree-matrix theorem) 



T spanning tree ofA eST 



(2.26) 



where the sum is over spanning trees on {1, ...n} oriented away from root i. 



Proof We give here a sketch of the Grassmann proof given in [39j. We can assume 
without loss of generahty that i = 1. For any matrix A we have: 



detA"= / []^rf?/'i(iV'i]?/'i?/'ie" 



4 = 1 



The trick is to use (2.25) to write 



J' 



hence 



det A 



by the Grassmann rules. We now expand to get 



n {l~A,,{^|J,-^P,)^|J,) 



detA 



11 



5 \^=(i,j)6g 



where Q is any subset of [n] x [n] , and we used the notation 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



Qg^ f d^d^ (^i^i) [ n [(^^ - I • (2-31) 

\{ij)6e / 

Then the theorem follows from the following 

Lemma 2.2 Qg = unless the graph Q is a tree directed away from 1 in which case 
fig = 1. 
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Proof Trivially, if belongs to Q, then the integrand of Qg contains a factor 
"ipi — "ipi = and therefore fig vanishes. 

But the crucial observation is that if there is a loop in Q then again Qg = 0. This 
is because then the integrand of flj^^n contains the factor 



1pT{k) - V'r(l) = (^r(fc) - ^r(fc-l)) H h (^r(2) " ^t{1))- 



(2.32) 



Inserting this telescoping expansion of the factor ipr{k) — V't{i) into the integrand of 
njF,7e, the latter breaks into a sum of {k — 1) products. For each of these products, 
there exists an a G Z/A;Z such that the factor {;ipr{a) — V'r(a-i)) appears twice : once 
with the + sign from the telescopic expansion of {ipT{k) ~ '4't{i)), and once more with 
a + (resp. — ) sign if (r(a;),r(a — 1)) (resp. (r(a — l),r(a))) belongs to JF. Again, 
the Grassmann rules entail that Qg = 0. □ 



To complete the proof of (2.26) every connected component of Q must contain 1, 
otherwise there is no way to saturate the dipi integration. 

This means that Q has to be a directed tree on {1, ...n}. It remains only to see 
that Q has to be directed away from 1, which is not too difficult. □ 

The interlude is over and we now turn to perturbative QFT and to the parametric 
representation of Feynman amplitudes. 



3 Parametric Representation of Feynman Ampli- 
tudes 

In this section we will give a brief introduction to the parametric representation of 
ordinary QFT on a commutative vector space M"^. We may take the example of 0^ 
bosonic theory but the formalism is completely general. 

3.1 Green and Schwinger functions in QFT 

In particle physics the most important quantity is the diffusion matrix S whose 
elements or cross sections can be measured in particle experiments. The S matrix 
can be expressed from the Green functions through the reduction formulas. Hence 
they contain all the relevant information for that QFT. 

These Green functions are time ordered vacuum expectation values of the fields 
0, which are operator-valued and act on the Fock space: 

GNizi, ...,zn) = (V'o, T[(j){zi)...(j){zN)]i'o)- (3.1) 
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Here ipo is the vacuum state and the T-product orders (j){zi)...(j){z^) according to 
increasing times. 

In the functional integral formalism the Green functions can be written as: 

Gr,{z,, zr,) - j^^jc^.^^.^^ • (3.2) 

Here C = Cq + Cint is the full Lagrangian of the theory. The Green functions 
continued to Euclidean points are called the Schwinger functions and are given by 
the Euclidean Feynman-Kac formula: 

N 

Sr,{z^,...,ZN) = Z-^ / J]0(z,-)e-^^(^("))""I^0, (3.3) 

i=i 

(3.4) 

For instance for the 0^ theory, Cint = ^\4>{x)'^ and we have 

^(0) = \d,md'm + \m<P{xf + j^<P{xf (3.5) 

where 



A is the (bare) coupling constant, which characterizes the strength of the in- 
teraction, the traditional factor 1/4! is inessential but slightly simplifies some 
computations. 

m is the (bare) mass, 

Z is the normalization factor, 

is an ill-defined "flat" product of Lebesgue measures '^'^(•^) Qa-oh. 
space time point. 



The coefficient of the Laplacian is set to 1 in (3.5) for simplicity. Although this 
coefficient actually in four dimensions flows through renormalization, it is possible 
to exchange this flow for a rescaling of the fleld 0. 

To progress towards mathematical respectability and to prepare for perturba- 
tion theory, we combine the e~ ^ '-''^'^^^''^'^^ Dcj) and the free normalization factor Zq = 
J e~ ■l^ ^'^'^^^^^'^^Dcj) into a normalized GauBian measure diJc{4>) which is well-deflned 
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on some subspace of the Schwartz space of distributions S'{R'^) [B]. The co- 
variance of this measure is the (free) translation invariant propagator C{x,y) = 
J (f){x)(j){y)dfic'{4') , which by shght abuse of notation we also write as C{x — y) and 
whose Fourier transform is 

In this way the Schwinger functions are rewritten as 

N 

Sm{zu...,zm) = Z-' / rr0(z,-)e-^«''^-*('^)ci/xc(0), (3.7) 

Jr^ i=i 

;-/Bd£.ntWx))<ix^^^^^^)_ (3.8) 

However this expression is still formal for two reasons; for typical fields the inter- 
action factor is not integrable over so that J^jCint{4') is ill-defined (infrared or 
thermodynamic problem) and in dimension more than 2 even when the interaction 
factor is restricted to a finite volume it is still ill-defined because for typical distri- 
butions 0, products such as (p'^i^x) are also ill-defined. This is the famous ultraviolet 
problem which requires renormalization (see |15]), but this problem is not addressed 
here, as we discuss solely the structure of the integrands in Feynman parametric rep- 
resentations, not the convergence of the integrals. The reader worried by ill-defined 
integrals in the rest of this paper for space-time dimension d larger than 2 should 
impose a ultraviolet regulator. This means he should replace replace C {p) by a better 
behaved Ck(p) such as 

p-K{p2+m2) ^-oo , , 



C.{p)= ,^ w 2^ , = / e-"(-+^W, (3.9) 



C,{x,y)= I e-°"^'-(^-^)'/^"4^. (3.10) 



so that 

poo 

n (rr ^ / ^-om2~(a:-j/)2/4Q!_ 

We now turn to perturbation theory in which the factor e~ ■I^r'' ^'"i^"^) is expanded as 
a power series. This solves the thermodynamic problem, at the cost of introducing 
another problem, the divergence of that perturbation expansion. This divergence 
which in the good cases can be tackled by constructive field theory |m HHl HZl HH] 
will not be treated in this paper. 
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Figure 3: A 0^ graph 



3.2 Perturbation theory, Feynman Graphs 

Wick theorem is nothing but the rule of pairing which computes the moments of a 
GauBian measure. It allows to integrate monomials of fields 

/ (P{xi)...(P{Xn)dflc{<P) = E n C(^V,^.e) (3.11) 

G e&G 

where the sum over G is over all contraction schemes (i.e. pairings of the fields) and 
C{xi^,Xj^,) is the propagator kernel joining the arguments of the two fields (t>{xi^) and 
4>{^je) pahed into the edge e by the contraction scheme G. 

It was Feynman's master stroke to represent each such contraction scheme by a 
particular graph in which edges represent pairs of contracted fields and vertices stand 
for the interaction. 

In the case of a 0"^ theory, remark that these interaction vertices have degree 4. 
Indeed the Schwinger functions after perturbative expansion are 

Sn{z^...zn) = / [ / Y['P\^v)dxv\^i.zi)..4{zN)dii{(p). (3.12) 

n=0 ^- J „=i 

The pairings of Wick's theorem therefore occur between n internal vertices each 
equipped with four fields and N external vertices or sources corresponding to the 
single fields 0(-2i), ... , 4'{zn). 

Schwinger function are therefore expressed as sums over Feynman graphs of as- 
sociated quantities or weights called the Feynman amplitudes. In this position space 
representation the Feynman graphs have both n internal vertices corresponding to 
the Cint factors, plus N external vertices of degree 1 corresponding to the fields 
(t>{z-i), (I){zn). In the case of the 0^ theory each internal vertex has degree 4. 
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The Feynman amplitudes are obtained by integrating over all positions of internal 
vertices the product of the propagator kernels for all the edges of the graphs 

Ag{zi,...,zn)^ Y[dxyYlC{xi^,Xj^), (3.13) 

V eeG 

where the product Yi^ runs over the internal vertices v. 

The quantities that are relevant to physical experiments are the connected Schwinger 
functions which can be written as: 



r,v(^i,...,^iv)= Yl ^-^^AiG)iz^,...,ZM), (3.14) 

4>^ connected graphs G with N{G)=N 



where S{G) is a combinatoric factor (symmetry factor). 

The momentum space representation corresponds to a Fourier transform to mo- 
menta variables called pi, ...,pn: 

riv(pi,...,Piv) = J dzi...dzNe'''^Pf'frN{zi,...,ZN), (3.15) 

where the factor 2 is convenient and we forget inessential normalization factors. 
This is a distribution, proportional to a global momentum conservation <^(X]7=iP/)- 
From now on we use an index / to label external momenta to remember that they 
are associated to corresponding graph-theoretic flags. Usually one factors out this 
distribution together with the external propagators, to obtain the expansion in terms 
of truncated amputated graphs: 



cf>^ truncated graphs G with N{G)=N 
N N 



-A)"(^) 



^(J2pf) n ;;2^^S(pi, -.pn). (3.16) 



/=1 



In this sum wc have to describe in more detail the truncated graphs G with N external 
flags. Such truncated graphs are connected, but they may contain bridges and self- 
loops. They no longer have external vertices of degree 1. Instead, they still have N 
external variables p/, no longer associated to edges but to flags {N in total), which 
decorate the former internal vertices. For instance for the (/)^ theory the degree of a 
truncated graph G is no longer 4 at each internal vertex. It is the total degree, that 
is the number of half-edges plus flags which remains 4 at every vertex. 
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Figure 4: A truncated 0"^ graph 



Ordinary Schwinger functions can be expressed as sums over partitions of the 
arguments of products of the corresponding truncated functions. We now give the 
exphcit form of the corresponding truncated amphtudes Aq{pi, ...,pn)- 



3.3 Parametric representation 

We shall first consider a fixed truncated oriented diagram G and compute the corre- 
sponding contribution or amplitude Aq as given by Feynman rules. 

We denote again by E and V the number of edges and vertices respectively, and 
by the number of fiags. Since G is connected its incidence matrix has rank V — 1. 

Now consider a Feynman graph G contributing to some truncated Schwinger 
function T'^{pi, ...,pn)- The usual way to take into account the global 6 function in 



(3.16) is to restrict to configurations such that ^fPf = 0. Extraction of this global 



delta function in (3.16) for the amplitude of a particular graph can be done provided 



we do not integrate the position of one of the vertices in (3.13), but rather fix it at 
an arbitrary point, eg the origin. From now on we suppose this vertex is v the one 
with last index. It provides a root in the graph G. However this standard procedure 
requires the non-canonical choice of that root vertex, and the final result does not 
depend on that choice. 

Another possibility is to modify the interaction A<^^(x) into Ae"''^ 4>^{x), in which 
case there is no longer global momentum conservation. One can compute modified 
amplitudes Bq(pi, ...pN^q) without factoring out the global S{J2'f=iPf) factor, so 
that 



r^(pi,...,P7v;g) = ^ 



^_A)«(G) 



SiG) 

4>'^ truncated graphs G with N{G)=N ^ ' 
^ 1 

Vi ^2 .^2 ^l^Vi.-.VN\q). (3.17) 
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The momentum conserving usual amplitudes are recovered when q ^ 0: 

N 

limSc(pi,...,p^;g) = (5( ...,p7v)- (3.18) 

a— >0 ' 



/=1 



This is the procedure we shall follow in subsection |3.4[ because it avoids the choice of 
a noncanonical root. But for the moment let us complete the standard presentation 
of Ag(pi, ...,P7V). 

The momentum representation of A^, forgetting from now on inessential factors 
of 271, is: 



v-i 



AI{p^,...,pn) = [f\d''KT^^^Y\6{ef,pf + eM. (3.19) 

in which we use the convention that repeated indices are summed, so that ef^pf+e^vke 
stands for the total momentum Ylif ^fvPf + J2e ^evke incoming at vertex v. 

To obtain the parametric representation we have first to rewrite the propagators 

as : 

1 

dae"'^^'+"''\ (3.20) 







J^2 _j_ ^2 

We obtain the momentum parametric representation 

„ E V-l 

AUpu -.Pn) = / l[daed%e-''^^''"+'^"^ J] 6{ef,pf + e^. (3.21) 

e=l v=l 

Fourier transforming the V — 1 Dirac distributions into oscillating integrals we 
obtain, up to some inessential global factors the phase-space parametric representa- 
tion 

„ E V-l 

AI{pi,...,pn) = / n [daee~''''"''d'^K] JJ ci'^a;^e-°'='='+2i{p/./.x„+fce.e„x„)^ ^3 22) 

e=l v=l 

where again keCeyXv means J2e=i Ylv=i ke^evXv etc, and the factor 2 is convenient. 

Finally integrating out the edge momenta whose dependence is GauBian leads to 
the X or direct space parametric representation: 



e=l D=l 
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Remember this amplitude is only defined on the submanifold = 0, because it is 
only there that the formula gives a result independent of the choice of the root not 



integrated out in (3.23) 



The parametric representation consists in integrating out fully the x or p variables 



in (3.21), (3.22) or (3.23). One obtains the parametric representation, which is an 



integral on a parameters only: 



E 



e=l 



-VG{p,a)/UG{a) 



(3.24) 



where Ug and Vg are called the first and second Symanzik's polynomials. 



Theorem 3.1 The first Symanzik polynomial Ug in (3.24) ^■^ ihe multivariate Tutte 



polynomial (2.10). On the submanifold pg = 0, the only one where it is unambigu- 



ously defined, the second polynomial Vg of (3.24) coincides with (2.13) and (2.15). 



We are going to give two proofs of this classic theorem of quantum field theory, 
one relying directly on contraction-deletion and on the phase-space representation 



(3.22) the other more standard and relying on the direct representation (3.23) and 



on the tree-matrix theorem. 

Indeed in order to compute the Symanzik's polynomials, let us remark first that 
the momentum representation mostly used in textbooks is not very convenient. To 



use (3.21) we should "solve" the 5 functions, that is rewrite each edge momentum 



in terms of independent momenta for cycles. In physics this is called a momentum 
routing. But such a momentum routing is linked to the choice of a particular spanning 
tree of G. The momenta of the edges not in this tree are kept as independent variables 
and the tree edges momenta are recursively computed in terms of those by progressing 
from the leaves of the tree towards the root which is the fixed vertex f„. This is not 
a canonical prescription, as it depends on the choice of the tree. 



The representations (3.22) or (3.23) are more convenient to integrate the space 



or momentum variables because the dependence in variables x and k is GauBian so 
that the result is a determinant to a certain power times a Gauf^ian in the external 



variables. In fact (3.22) is the best as we shall argue below. However there is still 



a small noncanonical choice, the one of the root. This is why we prefer to compute 
the regularized amphtudes 

„ s V 
Bj:{pi,...,pN;q) = / n [daee-''^"''d'^ke] JJ rf'^x^e""^^'-"^-! 

e=l v=l 

(3.25) 
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and to deduce the ordinary amplitudes from a limit g — 0. 

The last modification we perform is to attribute a different weight to each ver- 
tex regulator. This is more natural from the point of view of universal polynomials. 
So we define 

„ s V 

e=l v=l 

(3.26) 

These amplitudes are Gauf5ian in the external variables pf and no longer involve any 
noncanonical choice. We shall now compute their generalized Symanzik polynomials 
and deduce the ordinary Symanzik polynomials from these as leading terms when 
all g^'s are sent to 0. 



3.4 Generalized Symanzik Polynomials 



We consider the phase space representation (3.26). We have to perform a GauBian 
integral m. E + V variables (each of which is d-dimensional) . We consider these 
momentum and position variables as a single vector. We also forget the label for 
truncation as it is no longer needed in this section. The graph we consider may be 
connected or not. 

We introduce the condensed notations: 

BcipfAv) = llldaee-''^'^"d% IHd'^x^e-'^'''^''" (3.27) 

e V 

where Xq is a d{E + V + N) by d{E + V + N) square matrix, namely 

(tte ~ieev \ 
-itev qv -i<^fv (3.28) 
-lef, / 

where Oe and q^ are short notations for diagonal matrices tte^e.e' and q^Svy- Y is an 
+ + by 1 line, namely Y = [ke pj). 
We can further decompose Xq as 

^« - [X "f ) ■ i^-^) 

where Qg = ( '^^ ^^'^^ ) is a d[E + V) by d[E + V) square matrix and Rq is 

y ^Ceti Qv J 

the real rectangular N hy E + V matrix made of a dN by dE zero block and the dN 
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by dV "incidence flag" matrix e^^. The dimensional indices fi being quite trivial we 
no longer write them down from now on. 

Note P the line pf, hence the last part of the line Y. GauBian integrations can 
be performed explicitly and the result is a GauBian in external variables. Therefore 
up to inessential constants 

J „ det QA 



det Qq 
-\ 



p-V/U 

Wda^e-^^^^^'d^Kj-^ (3.30) 



for some polynomial Uq in a's and g's and a quadratic form in the p variable Vg 
with polynomial coefficients in a's and g's. 

Definition 3.1 The generalized Symanzik polynomials with harmonic regulators are 
the polynomials appearing in ( 3. 3(^ , namely 



l^ciae, qv) = det Qg, (3.31) 

VG(«e,g.,P/)/WG(«e,g.) = PRgQg^RgP'- (3-32) 

These polynomials can be computed explicitly: 
Theorem 3.2 

r e<i,T C 

VG{ae,qv,Pf) = ^Ylae^plYlqc', (3.34) 

where the sum over T runs over all forests of the graph, and the indices C and C 
means any connected component of that forest (including isolated vertices if any). 
The variables pc and qc are the natural sums associated to these connected compo- 
nents. 

In order to prove this theorem we introduce now the quasi- Grassmann represen- 



tations of Ug and Vg of Lemma 2.1 



Let's calculate first hence the determinant of Qg- Factoring out powers of i 
we get: 

detgc = det f"*^ (3.35) 
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which can be written as sum of a diagonal matrix D, with coefficients D^e = ote and 
^vv = Qv and of an antisymmetric matrix A with elements eev, that is, Q = D + A. 



By Lemma 2.1 



v,e 

Similarly V which is a minor related to the Qg matrix is given by a Grassmann 
integral but with sources 

v,e 

iXvUJv' + Xv'UJv), (3.37) 

where we have expanded ipvi'v' as ^[XvXv' + (jJv^v' + iiXv^v' + Xv'^v)] and canceled 
out the XvXv' + oJvOJv' term which must vanish by symmetry and the i factors. 

Now we can prove directly that these polynomials obey a deletion-contraction 
rule. 

Theorem 3.3 For any semi-regular edge e 

i^cic^e, Qv) = Oie UG-e{<^e, Qv) + i^G/e{oie, Qv) , (3.38) 
VG{(ye,qv,Pf) = aeVG-e{0ie,qv,Pf) + VG/e{(ye,qv,Pf)- (3.39) 

Moreover we have the terminal form evaluation 

UG^ae, qv) = n «e n {"iAQ) 

e V 

VG{ae,qv,Pf) = Wae^plWqv (3.41) 

e V v'^v 

for G solely made of self-loops attached to isolated vertices. 

Proof If G is not a terminal form we can pick up any semi-regular edge e connecting 
vertices Vi and V2 with e^,^ = +1, 6^,3 = —1. We expand 

g-aeXe<.. _ I + (3,42) 
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For the first term, since we must saturate the Xe and u^, integrations, we must keep the 
Xe(Xfi ~Xv2) term in 6^"-*^='^'=''^" and the similar uj term, hence we get a contribution 

Performing the trivial triangular change of variables with unit Jacobian: 

Xi-i = Xi-i - Xi>2> Xv = Xv for v^vi, (3.44) 

and the same change for the u variables we see that the effect of the {xvi —Xv2){^vi — 
LOv^) term is simply to change the vi label into V2 and to destroy the edge e and the 
vertex vi. This is exactly the contraction rule, so det QG,e,i = detQc/e- The second 
term det QG,e,2 with the aeOUeXe factor is even easier. We must simply put to all 
terms involving the e label, hence trivially detQG,e,2 = OedetQc-e- Remark that 
during the contraction steps the weight factor qviXvi^vi is just changed into qviXv2^v2- 
That's why we get the new weight + qv2 for the new vertex V2 which represent 
the collapse of former vertices vi and V2- 

Note that the source terms in V do not involve Xe and Ue variables. Therefore 
the argument goes through exactly in the same way for the second polynomials. The 
only remark to make is that like weights, flag momenta follow contraction moves. 

The evaluation on terminal forms is easy. For a graph with only vertices and self- 
loops the matrix Qg is diagonal, because eev is always 0. Hence Ug is the product of 
the diagonal elements He "^e lit; 1v The second polynomial can be analyzed through 



the Grassmann representation, but it is simpler to use directly (3.32) and the fact 



that Qg is diagonal to get (3.41). This completes the proof of Theorem 3.3, hence 
also of Theorem 13.21 □ 

We turn now to the limit of small regulators to show how for a connected graph 
G the ordinary amplitude pf)AG and the ordinary polynomials Ug and Vg 

emerge out of the leading terms of the regularized amplitude Bg and the generalized 
polynomials Ug and Vg- 

When all g's are sent to zero there is no constant term in Ug but a constant term 
in Vg- Up to second order in the q variables we have: 

UG{ae,qv) =qGj2Il^^ + '^^^^'>^ (^.45) 
r e^r 

VG{ae,qv,Pf) = pI^YI^^ + ^{ph,qG2 +Pg2^g^) n (^e + 0{q^), (3.46) 
r e^r r2 e^r2 
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where the sum over T runs over trees and the sum over 7^ runs over two trees 
separating the graph into two connected components Gi and G2- Hence we find 



-v/u 



-Pg/IG g-ET2(PGi''G2+PG2 9Gi)ne5^r2°<=/9G Er Oe^^T "<=+P^O(l)+0(g) 



d/2 

<1g 



(3.47) 



Up to inessential normahzation factors the first term tends to S{pg) and the second 
one tends to e~'^/^ /U'^/'^ if we use the fact that 5{pg) f [pg] = '^(pg)/(0), that is if we 
use the delta distribution to cancel the PqO(1) term and to simplify {PG1QG2+PG2IG1) 
into QcPgi = 1gPg2- This proves (3.18). 

The Ug and Vg polynomials are in fact easy to recover simply from the Ug 
polynomial alone: 



Theorem 3.4 For any connected G and any vertex v 

d 



UGiae 



dqv 



UGiae, Qv) 



\q^,=0 W- 



On the submanifold pg = we further have 



VG{ae,Pf) 



Pv 



dqvdq. 



-UG{ae,qv) 



1 5„" =0 Vi;" ■ 



Proof It is an easy consequence of Theorem 3T2| 

We can also prove an analog of Proposition 
on the submanifold pg = 0- 



(3.48) 



(3.49) 



□ 



2.1 



between Vg and IAg{vv') but only 



3.5 Relation to discrete Schrodinger Operator 

As an aside, it is worthwhile to notice that there is a relation with discrete Schrodinger 
operators on graphs f53|. Recall that given a graph G = (y,E), the discrete Lapla- 
cian is defined as follows. We first introduce the 0- forms flo{G) = MX as the real 
functions on the set of vertices and 1-forms ^lo{G) = as functions on the edges. 
Then, the discrete differential d : ^lo{G) ^i{G) is defined as 

d^(e) = ^ee^,^„, (3.50) 

V 

where we recall the convention that for a self-loop eev = and an arbitrary orientation 
is chosen on the edges. Next, given strictly positive weights Pe associated to the edges, 
we define d* : (]i(G) ^ fio(G') by 

d*cj){v)=J2l3eeev(f>e. (3.51) 
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Note that d* is precisely the adjoint of d for the scalar product on defined by 
the weights Pe and the Euclidean one on MX . Accordingly, the 0-form Laplacian 
A : no{G) no{G) is 

A = d*d, (3.52) 
or, in terms of its action on functions ip G MX , 

(3.53) 

e,v 

Note that there is exactly one zero mode per connected component, as follows from 
the equivalence between A?/^ = and Aifj = Finally, the weights associated 
to the vertice^ define a function V from the vertices to M acting multiplicatively 
on VLq{G) so that we define the discrete Schrodinger operator (Hamiltonian in the 
quantum mechanics language) on the graph by 

H = -A + V. (3.54) 

Turning back to the parametric representation, if we perform the GauBian integration 
over the momenta we are left with 



D/2 



(ai ■ ■ -ae)^/' 



J Yl (ix^e"^"'""^"''^"'"'^'''"'"^*^'''^"'^", (3.55) 



with weights /?e = In particular, the first Symanzik polynomial with regulators 
g„ is expressed in terms of the determinant of H, 

WG(«,g) = (rie^e) detH= (n,ae) j U.d^J^ve-^^'^'^^''^y^^', (3.56) 

with ip^,ipv Grassmann variables. By the same token, the ratio appearing in the 
Feynman amplitude is expressed in terms of its inverse G (Green's function in the 
quantum mechanics language). 



VG{a,q,p) 



J^^G^yPi, (3.57) 



where the Green's function can also be expressed using Grassmann integrals. As a 
byproduct, it turns out that it can also be computed by contraction/deletion. 



^Strictly speaking, the latter are associated to the flags and g„ is the sum the weights of the 
flags attached to v. 
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3.6 Categorified Polynomials 



We have up to now considered two seemingly unrelated graph polynomials obeying 
contraction/deletion rules, the multivariate Tutte polynomial ZciPe, q) and Uaic^e, Qi), 
from which the Symanzik polynomials can be recovered by various truncations. 
Therefore, it is natural to wonder wether there is a single graph polynomial, obey- 
ing contraction/deletion rules too, from which both ZG{f3f,,q) and UciacQi) can be 
recovered. In this subsection for simplicity we shall consider only the first Symanzik 
polynomial, and the flags considered in this subsection no longer bear external mo- 
menta, but an abstract index. 

Such a polynomial is an invariant of graphs with flags, i.e. labeled half-edges 
attached to the vertices. In order to make the contraction possible, it is necessary 
to allow each vertex to have several flags, all carrying distinct labels. The requested 
polynomial, WciPe, Qi) depends on edge variables Pe as well as on independent vari- 
ables qi for each non empty subset / of the set of labels of the flags, with the proviso 
that, for each vertex, the subsets / contain all the flags attached to the vertex or 
none of them. Thus, for a diagram with V vertices carrying flags there are 2^—1 
variables qj. 

Definition 3.2 For a graph G with flags, WaiPe-i Qi) is defined by the expansion 

woiPe: = E ( n n ' ^^-^^^ 

connected components 

where In are the sets of flags attached to the vertices of the connected component Cn 
of the spanning graph iV,A). 

For example, for the bubble graph on two vertices with two edges between these 
vertices and flags 1, 2 attached to one of vertex, and flag 3 to the other one, we have 

WaiPe, qi) = (/9i/92 + /9i + /92)gi23 + gi2g3. (3.59) 

Since the variables qi are defined using the flags, the contraction/deletion rule for 
yVciPe: Qi) requires us to properly dcflne how the flags follow the contraction/deletion 
rule for any edge of G — e and G/e. Because the vertices and the flags of G — e are 
left unchanged, the same variables qi appear in G and G — e. For G/e, we restrict 
the qi to those associated with subsets that contain either all the flags attached to 
the two vertices merged by the contraction of e, either none of them. This is best 
formulated using flags: the new vertex simply carries the flags of the two vertices 
that have been merged. Then, the contraction/deletion identity simply follows from 
grouping the terms in WciPe, Qi) that contain /3e and those that do not. 
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Proposition 3.1 The 'polynomial WG(/?e, (li) obeys the contraction /deletion rule for 
any edge 

WciPe, qi) = Pen^G/eiPe'^e, qilc/e) + Wc-eiPe'^e, Ql). (3.60) 

The multivariate Tutte polynomial is easily recovered by setting qi = q for any 

^G(/5e,g) = >VG(/3e,g/ = g). (3.61) 

In this case, all the information about the flags is erased and so that the latter may 
be omitted. To recover UG{ae,qi), it is convenient to introduce as an intermediate 
step the polynomial 

TG(«e, ?o = 5^ n n ( E ^0 ' ^^-^^^ 

AcE e^E Cn i^In 

where as before /„ are the flags included in the connected component C„ of the 
spanning graph {y,A). By its very definition, Tcioie^qi) is related to 'WciPeiqi) by 
setting qi = Y^iaiq^ 

^G{ae,qt) = (n^^) ^G{Pe = l/(ye,qi = Y.iei^i)- (3-63) 

e 

Then, the polynomial Ucideyqi) is obtained from Tcidcqi) by keeping only the 
highest degree terms in the a/s for each term in Y[c„ Sie/n Indeed, Uciae, qi) 
is obtained from TG{ae,qi) by truncating its expansion to those subsets A C E 
that are spanning forests, i.e. that obey = |A| — V + k{A). Since the number of 
connected components k{A) is fixed by the global degree in the g^'s, the forests are 
obtained with 1^41 minimal, so that the global degree in the a^s must be maximal. 
Note that a truncation to the spanning forests may also be performed at the level of 
the multivariate Tutte polynomial by restricting, at fixed degree in q, to the terms 
of minimal degree in the P/s. This yields an expansion over spanning forests |1] (see 
also [12]). 

FG{Pe,q)= Yl (H/^O^'^^^- ^^-^^^ 

ACE eeA 

spanning forest 

This, as well as the relation to the Symanzik polynomial, is conveniently summarized 
by the following diagram. 
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Proposition 3.2 The previous polynomials may be obtained from yV{ae,qi) by the 
following series of substitutions and truncations, 



highest order 
\in the Op 



G{ae,qi) ^UG{ae,qi) (3.65) 

multiplication by Yle'^e ^^^^^ term in 2_^i Qi 





Qt — (7^\ , , , term inq 

' lowest order ^ , , , ^„ 

^ ■ ^ multiphcation by l^ctp 

ZG{Pe,q)'^Fa{Pe,q) 
where ae = l/Pe- 

Alternatively, the polynomial WcicMe, qi) can be seen as an extension of the poly- 
nomial Wciiaiy) introduced by Noble and Welsh in |13]. 

Definition 3.3 For a graph with weights uj^ G N* assigned to the vertices, the W 
polynomial is defined as 

Wc{ia.y) = Y. (y-l)'-"-^^^) n (3.66) 
Ace Ci,...,Ch^A) 

connected components of{V.A) 

with a„ = ^^gc ^'^''^ of the weights of the vertices in the connected component 

C 

This polynomial also obeys the contraction/deletion rule if we add the weights 
of the two vertices that are merged after the contraction of an edge. Alternatively, 
weights may be assigned to flags, with the convention that the weight of a vertex 
is the sum of the weights of the flags attached to it. Then, W^^a^u) is naturally 
extended to diagrams with flags and results from a simple substitution in Wciia, y)- 

Proposition 3.3 For a graph with weights uoi G N* assigned to the flags, 

WaiUy) = (y - lp''^WG{Pe=y-l,qi = iy -l)^a,), (3.67) 

with aj = ^he sum of the weights of the flags in I . 

The polynomial Wciia^y) only encodes the sum of the weights of the flags in 
each connected component and erases information about their labels. In particular, 
if we weight each flag by cjj = 1, then the expansion of W only counts the number 
flags per component whereas that of WG(/?e, qi) keeps track of the associated set of 
labels. In a more sophisticated language, the latter may be considered as the simplest 
categorification of the former: integers, understood as finite sets up to isomorphisms, 
have been replaced by the category of finite sets. 
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3.7 Symanzik Polynomials through the tree matrix theorem 
in x-space 



In this section we provide a sketch of a more standard proof of Theorem 3.1 through 
the X space representation and the tree matrix theorem. The reason we include it 
here is for completeness and because we have not been able to find it in the existing 
literature, in which the same computation is usually performed through the Binet- 
Cauchy theorem. 

The V X V matrix Qg{c() analog in this case of (3.29) is defined as 

[QG{a)]v,v' = y^eev—eev'- (3.68) 

e 

It has vanishing sum over lines (or columns): 

'^[QG{oi)]vy = XlXl ^^e— W = 0. (3.69) 

V' v' e 

Therefore by the tree matrix theorem the determinant of the {V ~ 1) x (V — 1) 
matrix Qg{(^) defined as its principal minor with the line and column for the root 
vertex number V deleted is: 

Ada) = det[QG(«)] =$^11- (^.70) 

r eeT 

where the sum is over all trees of G. Since every tree of G has V — 1 edges, Aq 
is clearly a homogenous polynomial in the a~^. For a > 0, A is positive. The 
remaining {V — 1) vectors z may then be integrated over and the result is 

MP) = /'°°n(^°^^-°-') "''i''''°"^A''r''!;'r^ - 

Jo V [ai...aEAG{a)\'^/^ 

This formula expresses Ag{p) as a function of the invariant scalar product of 
external momenta p„ ■ The denominator 

Ucia) = ai-.-aEAcia) = ^Y\_^e (3.72) 

r e^r 



is a homogenous polynomial of degree V—1. This gives an alternative proof of (3.24 ). 
The second Symanzik polynomial can also be obtained through this method and the 
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corresponding computation is left to the reader. Of course harmonic regulators can 
also be included if one wants to avoid the noncanonical choice of a root, but the 
Pfaffian structure of the phase space representation is lost. Also this x-space method 
does not generalize easily to noncommutative field theory to which we now turn our 
attention. 

4 Bollobas-Riordan Polynomials 
4.1 Ribbon graphs 

A ribbon graph G = (V, E) is an orientable surface with boundary represented as 
the union of V closed disks, also called vertices, and E ribbons also called edges, 
such that: 

• the disks and the ribbons intersect in disjoint line segments, 

• each such line segment lies on the boundary of precisely one disk and one 
ribbon, 

• every ribbon contains two such line segments. 

So one can think of a ribbon graph as consisting of disks (vertices) attached to each 
other by thin stripes (edges) glued to their boundaries (see Figures [5||6]) . For any 
such ribbon graph G there is an underlying ordinary graph G obtained by collapsing 
the disks to points and the ribbons to edges. 

Two ribbon graphs are isomorphic if there is a homeomorphism from one to the 
other mapping vertices to vertices and edges to edges. A ribbon graph is a graph 
with a fixed cyclic ordering of the incident half-edges at each of its vertices. 

A face of a ribbon graph is a connected component of its boundary as a surface. 
If we glue a disk along the boundary of each face we obtain a closed Riemann surface 
whose genus is also called the genus of the graph. The ribbon graph is called planar 
if that Riemann surface has genus zero. 

Generalized ribbon graphs that can also incorporate Moebius strips and corre- 
spond to nonorientable surface can be defined but will not be considered in this 
paper. 

There is a duality on ribbon graphs which preserves the genus but exchanges 
faces and vertices, keeping the number of edges fixed. It simply considers the disks 
glued along faces as the vertices of a dual graph and changes the ends of each ribbon 
into borders of the dual ribbon. 
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Extended categories of ribbon graphs with flags can be defined. Flags can be 
represented as ribbons bordered by dotted lines to distinguish them from ordinary 
edges (see Figures [5||6]) . Beware that the cyclic ordering of flags and half-edges at 
each vertex is very important and must be respected under isomorphisms. The genus 
of an extended graph is defined as the genus of the graph obtained by removing the 
flags and closing the corresponding segments on their vertices. The number of broken 
faces is the number of faces which do contain at least one flag. It is an important 
notion in noncommutative field theory. 

We define for any ribbon graph 

• V{G) as the number of vertices; 

• E{G), the number of edges, 

• k{G), the number of connected components, 

• r{G) = V{G) - k{G), the rank of G, 

• n{G) = E{G) - r(G), the nullity of G, 

• bc{G) = F{G), the number of components of the boundary of 

• g{G) = k — (y — E + be) /2 is the genus of the graph, 

• f{G) the number of flags of the graph. 

A graph with a single vertex hence with = 1 is called a rosette. 

A subgraph if of a ribbon graph G is a subset of the edges of G. 

The Bollobas-Riordan polynomial, which is a generalization of Tutte polynomial, 
is a algebraic polynomial that is used to incorporate new topological information 
specific to ribbon graphs, such as the genus and the number of "broken" or "external" 
faces. It is a polynomial invariant of the ribbon graph. 

4.2 Bollobas-Riordan Polynomial 

Definition 4.1 (Global definition) The Bollobas-Riordan polynomial is defined 
by: 

= R^(x, y,z)=J2i^- (4.1) 

HcG 

•^This is the number of faces of G when G is connected. 
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Figure 5: A planar ribbon graph with V — E — 1. be — 2 and two flags. 




Figure 7: Contraction of the single self-loop Gi. 



The relation to the Tutte polynomial for the underlying graph G is Rg{x — l,y — 
1,1) = TQ{x,y). Remark also that if G is planar we have Rg{x — l,y — 1, z) = 

When if is a spanning graph of G, we have k{H) — k{G)=r{G) — r{H). So we 
can rewrite the R polynomial as: 

Rg = {x- 1)-'=(^) M{H), (4.2) 

HCG 

where 

MiH) = ix- l)HH)yniH)^kiH)-bc(H)+n(H) (4 3) 

SO that M{H) depends only on H but not on G. 



4.3 Deletion/contraction 

The deletion and contraction of edges in a ribbon graph are defined quite naturally: 
the deletion removes the edge and closes the two scars at its end; the contraction of a 
semi-regular edge creates a new disk out of the two disks at both ends of the ribbon 
with a new boundary which is the union of the boundaries of the two disks and of 



the ribbon (see Figure 10). An interesting property is that deletion and contraction 
of edges are exchanged in the dual graph. 

The deletion of a self-loop is standard. However the natural contraction of a 
self-loop creates a surface with a new border. Iterating, we may get surfaces of 
arbitrary genus with an arbitrary number of disks removed, a category also called 
disk-punctured surfaces. The ribbons can now join any puncture to any other. For 
instance the contraction of the self- loop on the graph Gi of Figure [7] leads to a cylinder 
ie to a single vertex which is a sphere with two disks removed. The contraction of 
the two self-loops in graph G2 of Figure [8] corresponds to the cylinder with a ribbon 
gluing the two ends, hence to a torus. 
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Figure 8: Contraction of the two self loops non-planar G2- 



Deletion and contraction defined in this extended category of graphs can be it- 
erated until the graph has no longer any edge, ie is a collection of disk-punctured 
Riemann surfaces. These punctured Riemann surfaces are very natural objects both 
in the context of string theory and in NCQFT. However we do not consider them in 
this paper. 

In this paper we remain in the category of ordinary ribbon graphs with disk-like 
vertices. The contraction/deletion of semi-regular edges leads to rosettes as terminal 
forms. To treat them we introduce the notion of double contraction on nice crossings. 
Nice crossings were introduced in [20] : 

Definition 4.2 A nice crossing pair of edges in a rosette is a pair of crossing edges 
ci and 62 which are adjacent on the cycle of the rosette. Adjacency means that one 
end of ei is consecutive with an end of 62 (see Figure^. 

It is proved in [20] that any rosette TZ of genus g > contains at least one nice 
crossing. 

The double contraction of such a nice crossing pair consists in deleting ei and 62 
and interchanging the half-edges encompassed by ei with the ones encompassed by 
62, see Figure [9] This double contraction was defined in [20] under the name of "3rd 
Filk move". It decreases the genus by one and the number of edges by 2. 

In the next section iterating this double contraction until we reach planarity 
allows us to compute the U* Symanzik polynomial by remaining in the category of 
ordinary ribbon graphs. 

Theorem 4.1 (Bollobas-Riordan polynomial, contraction/deletion) 

Rg = Rc/e + Rc-e (4.4) 
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Figure 9: When deleting the two edges of a nice pair crossing on some contracted 
vertex, one also needs to interchange the half-edges encompassed by the first edges 
with those encompassed by the second one. Beware that the horizontal line in this 
picture is a part of the rosette cycle. 




Figure 10: The contraction-deletion for a ribbon graph. 
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for every ribbon graph G and any regular edge e of G and 



Rg = xRc/e (4.5) 

for every bridge of G. 

Therefore the R polynomial satisfy contraction-deletion relations as the Tutte 
polynomial. However to complete its definition we also need to define the R polyno- 



mial for single vertex graphs, namely the rosettes, which can be read off from (4.3). 
For such a rosette TZ, k{TZ) = V{TZ) = k{H) = V{H) = 1, so that the R polynomial 
does not depend on x and 

Rn{y,z)=J2y^^^h'^^^\ (4.6) 
Hen 

For 2; = 1 we recover Rfi{y — 1,1) = y^^'^K 

4.4 The multivariate Bollobas-Riordan polynomial 

Like in the case of Tutte polynomial, we can generalize the Bollobas-Riordan poly- 
nomial to a multivariate case. As before,we associate to each edge e a variable (3e- 

Definition 4.3 The multivariate Bollobas-Riordan polynomial of a ribbon graph 



analog of the multivariate polynomial |^.^[ ) is: 

ZG{x,{Pe},z) = J2 x"'''\llPe)z''^''\ (4.7) 



HCG eeH 



It obeys again a deletion/contraction relation similar to Theorem (4.1) for any 
semi-regular edge. 



5 Translation-invariant NCQFT 
5.1 Motivation 

Noncommutative quantum field theory, hereafter called NCQFT, has a long story. 
Schrodinger, Heisenberg [S^ and Yang [55] tried to extend the noncommutativity 
of phase space to ordinary space. Building on their ideas Snyder [56] formulated 
quantum field theory on such noncommutative space in the hope that it might behave 
better than ordinary QFT in the ultraviolet regime. 
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Right from the start another motivation to study noncommutative quantum field 
theory came from the study of particles in strong magnetic fields. It was early 
recognized that non zero commutators occur for the coordinates of the centers of 
motion of such quantum particles, so that noncommutative geometry of the Moyal 
type should be the proper setting for many body quantum physics in strong external 
field. This includes in condensed matter the quantum Hall effect (see the contribution 
of Polychronakos in [SZ]), or other strong field situations. 

An other motivation comes from particle physics. After initial work by Dubois- 
Violette, Kerner and Madore, Connes, Lott, Chamseddine and others have force- 
fully advocated that the classical Lagrangian of the current standard model arises 
naturally on a simple noncommutative geometry. For a review see Alain Connes's 
contribution in [ST] and references therein. 

Still an other motivation came from the search of new regularizations of non- 
Abelian gauge theories that may throw light on their difficult mathematical structure. 
After 't Hooft proposed the large N limit of matrix theory, in which planar graphs 
dominate, as relevant to the subject [SB], the Eguchi-Kawai model was an important 
attempt for an explicit solution. These ideas have been revived in connection with 
the ultraviolet behavior of NCQFT on the Moyal- Weyl geometry, which also leads 
to the domination of planar graphs. Seiberg and Witten proposed in [59] a mapping 
between ordinary and noncommutative gauge fields which does not preserve the 
gauge groups but preserve the gauge equivalent classes. 

The interest for non commutative geometry also stems from string theory. Open 
string field theory may be recast as a problem of noncommutative multiplication of 
string states [60]. It was realized in the late 90's that NCQFT is an effective theory 
of strings [6T]. Roughly this is because in addition to the symmetric tensor g^^iy the 
spectrum of the closed string also contains an antisymmetric tensor B^,^. There is no 
reason for this antisymmetric tensor not to freeze at some lower scale into a classical 
field, inducing an effective non commutative geometry of the Moyal type. There 
might therefore be some intermediate regime between QFT and string theory where 
NCQFT is the relevant formalism. The ribbon graphs of NCQFT may be interpreted 
either as "thicker particle world-lines" or as "simplified open strings world-sheets" 
in which only the ends of strings appear but not yet their internal oscillations. 

5.2 Scalar models on the Moyal space 

The noncommutative Moyal space is defined in even dimension d by 

[x^a;^]. = ^e^^ (5.1) 



41 



where is an antisymmetric d/2 by d/2 block-diagonal matrix with blocks: 



e 

-9 



and we have denoted by -k the Moyal-Weyl product 



(5.2) 



{f^9){x) = I 0^,d'yf{x+le-k)g{x + y)e"'-y. (5.3) 

Note that in the limit 6^0 this product becomes the ordinary commutative product 
of functions. 

5.2.1 The "naive" model 

The simplest field theory on this space consists in replacing the ordinary commutative 
local product of fields by the Moyal-Weyl product 

5[0]= /"c/^x(^9^0^9'^0+^/iV* 0+^0*0*0*0)- (5.4) 



In momentum space the action (5.4) writes 



S[4>] = j d''p{^p^<Ppy + + V(0, 9)). (5.5) 

where V{(f), 9) is the corresponding potential. 

An important consequence of the use of the non-local product -k is that the inter- 
action part no longer preserves the invariance under permutation of external fields. 
This invariance is restricted to cyclic permutations. Furthermore, there exists a ba- 
sis - the matrix base - of the Moyal algebra where the Moyal-Weyl product takes 
the form of an ordinary (infinite) matrix product. For these reasons the associated 
Feynman graphs are ribbon graphs, that is propagators should be drawn as ribbons. 

In j62] several contractions on such a Feynman graph were defined. In particular 



the "first Filk move" is the contraction introduced in subsection 4.3 Repeating this 



operation for the V — 1 edges of a spanning tree, one obtains a rosette (see Figure 



nj). 

Note that the number of faces or the genus of the graph does not change under 
contraction. There is no crossing between edges for a planar rosette. The example 



of Figure 11 corresponds thus to a non-planar graph (one has crossings between the 



edges ki and k2). This pair is called a nice crossing pair. 
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Figure 11: An example of a rosette with two flags. The crossings of edges ki and k2 
indicate the non trivial genus (here g = 1). 

The notions expressed in the previous section (namely the Green and Schwinger 
functions or the perturbation theory concepts) remain the same as in QFT. Usual 
Feynman graphs are simply replaced by ribbon Feynman graphs. 



Recall that this "naive model" (5.4) is not renormalizable in ci = 4. This is due to 
a new type of non-local divergence at the level of the 2— point function - the UV/IR 
mixing |19j . 



5.2.2 A translation-invariant renormalizable scalar model 

In order to restore renormalizability at d = 4, the propagator can be modified in the 
following way [63] 

Sgmrt[4'] = J d'^P {\pf.^P^(t> + ^rri^H + +^0*0*0*0), (5.6) 

where a is some dimensionless parameter which is taken in the interval < a < 



The corresponding propagator writes in momentum space 



Cgmrt - o , ,,2 , ^ ■ (5- 

In [63] , this model was proved to be renormalizable at any order in perturbation the- 
ory. Furthermore, its renormalization group fiows [M] were calculated; a mechanism 
for taking the commutative limit has been proposed [65] (for a review on all these 
developments, see [UB]). 
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5.3 The NC Parametric representation 

In this subsection we present the implementation of the parametric representations 
for the noncommutative scalar models introduced in the previous subsection. 

To keep track of the cyclic ordering at the vertex it is convenient to detail the 
incidence matrix Eev into a more precise incidence tensor where i = 1 , . . . , 4 indexes 
the four corners of the Moyal vertex. As before it is 1 if the edge e starts at corner 
i of vertex w, -1 if it exits at that corner, and otherwise. 



To implement the parametric representation we follow subsection 3^ The prop- 
agator remains the same as in QFT, but the contribution of a vertex v now corre- 
sponds to a Moyal kernel. In momentum space it writes using again summation over 
repeated indices 



i=l 



By kiQkj we denote k'^Q^^k^-. The 5— function appearing in the vertex contribution 



(5.8) is nothing but the usual momentum conservation. It can be written as an 
integral over a new variable x^,, called hyperposition. One associates such a variable 
to any Moyal vertex, even though this vertex is non-local: 

Sii:^:^) = J (f^e-^Et^<»^^) = J i^e^-(S^-<^^^). (5.9) 

where a is a (i/2 by (i/2 block-diagonal matrix with blocks: 

. = (° ■ (5.10) 



Note that to pass from the first to the second line in (5.9), the change of variables 
ix[ = x^a has Jacobian 1. 

5.3.1 The "naive" model 

Putting now together the contributions of all the internal momenta and vertices, one 
has the following parametric representation: 

E 

AUpu ...,Pn) = K^ n dXda^e-'^^^'^'^"^'^ (5.11) 

e,e'=l 



y-1 



v=l 
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where we have denoted by some inessential normahzation constant. Furhermore 
note that in the integrand above we have denoted, to simphfy the notations, by ke 
or ke' momenta which can be both internal or external. 



5.3.2 The translation-invariant model 



The parametric representation of the model (5.6) was analyzed in [67]. This rep- 
resentation is intimately connected to the one of the model (5.4) (see the previous 

as 

(5.13) 



subsubsection) for the following reason. One can rewrite the propagator (5.7) as 

1 _ 1 
A + B ~ ~A 



-B- 



for 



A 



p + m , 



A A + B 



B 



(5.14) 



Thus, the propagator (5.7) writes 

1 



1 



GMRT 



p2 _|_ ^2 



p2 _|_ ^2 Q2p2(jp, _|_ j^2-) _|_ q' 

1 a 



p2 _|_ j^2 p2 _|_ j^2 Q2(jp, _|_ ^2^(j92 _|_ j^2j 



(5.14) 



where —m\ and —m\ are the roots of the denominator of the second term in the 
LHS (considered as a second order equation in p^, namely -^^"^^^fc^^^™'^-^^^" < g. 
Note that the form (5.15) allows us already to write an integral representation of 
the propagator C(p, m, 9). Nevertheless, for the second term one would need a triple 
integration over some set of Schwinger parameters: 



GMRT 



rf«e-"(^'+'"'), (5.15) 



a 

Instead of that one can use the following formula: 
11 1.1 



1 



rrin 



+ 1712 



This allows to write the propagator (5.15) as 

1 a 



1 



GMRT 



p2 _|_ ^2 ^2^^2 _ j^2^ p2 _|_ j^2 p2 _|_ ^2 p2 _|_ ^2 - 
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(5.15) 



(5.16) 



This form finally allows to write down the following integral representation: 

nOO POO POO 

Cgmrt = / rfae-"(p'+"^') - — ^ ^ / / darfaie'^^+^^^P'-"™'-"^'-? 

Jo ^^{mi -mf) Jq Jo 

POO POO 

+ -T^TT—o ^ / / rfac/a2e-("+"^)p'e-"'^'e-"2'"2. (5.16) 

d^{mi - m{) Jo Jo 

Let us also remark that the noncommutative propagator Cgmrt is bounded by 
the "usual" commutative propagator C{p,m) 

Cgmrt <Cip,m). (5.17) 



Using now (5.15), the parametric representation of the model (5.6) is thus a sum 
of 2^ terms coming from the development of the E internal propagators. Each of 
these terms has the same form of the one of polynomials in the previous subsection. 
The only differences comes from 

• the proper substitution of the set of Schwinger a parameters 

• the mass part. 
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One has 



+(--)^ 

ii<j2,iij2=i s 

1 



I<j2,ilj2 = l j7^jl,i2,i=l 



+ . . . + 



t7 J ~_T uyio:, ■ + ' + ' 112 



5.4 Deletion/contraction for the NC Symanzik polynomials 

In this subsection we give some results relating the Bollobas-Riordan polynomial 
and the parametric representations of the noncommutative scalar models introduced 
here. 



5.4.1 The "naive" model 



As in the commutative case, we have to perform a GauBian integral in a d{E + V — 1) 
dimensional space. Consider a ribbon graph G with a root v. 



We introduce the condensed notations analog to (3.27)-(3.29) 



Ag{p) 



Wdae 



-aem 



d'^xd'^pe 



-YXY^ 



(5.11) 
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where 



Y= {k 



e Pe 



X 



X 



-iR M 



(5.12) 



Q is an d{E + V — l)-dimensional square matrix. We have denoted by Pe the external 
momenta and by Xjj the hyperposition associated to the root vertex v. The matrix 
i? is a d{N + 1) X d{E + — 1) dimensional matrix and M is a d{N + 1) dimensional 
square matrix representing the Moyal couplings between the external momenta and 
the root vertex. 

GauBian integration gives, up to inessential constants: 



Y\dae 



det 



(5.13) 



where P is a line matrix regrouping the external momenta (and the hyperposition 
associated to the root vertex). 

The determinant of the matrix Q defines therefore the first Symanzik NC-polynomial 
U* and the product of the matrices R and inverse of Q defines the quotient of the 
second Symanzik polynomial V* by U* where the star recalls the Moyal product used 
to define this NCQFT. 

Let us calculate first the determinant of Q. One has 



(5.14) 



where D is a diagonal matrix with coefficients D^e = oie-, for e = 1, . . . , £" and D„y = 
for the rest, f = l,...,V^ — 1. A is an antisymmetric matrix. In |20j it was noted 
that, for such a matrix 



detQ = det{D + AY. 
This implies, as in the commutative case, that 

U"" = det{D + A). 
Factoring out powers of i one has 



(5.15) 



(5.16) 



det{D + A) = det 



/ "1 

-fl2 



fl2 

Oi2 



(5.16) 
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The difference witli tlie commutative case comes from tlie non-trivial antisymmetric 
coupling between the E edges variables. It corresponds to an E dimensional square 
matrix F with matrix elements 

n 4 

fee' = E E ^(^'^>eA' Ve < c', 6, c' = 1, . . . , (5.17) 

v=l ij=l 

where u is an antisymmetric matrix such that uj{i,j) = 1 if i < j. This matrix takes 
into account the antisymmetric character of G in k^j_Q^'^py. 



Using again Lemma 2.1 



dei{D + A) = / Y\_dujidxidujedxe 

i,e 

Note that the last term above represents the difference with the commutative case. 



We have the exact analog of Theorem 3.3 to prove a deletion-contraction rule. 



Theorem 5.1 For any semi-regular edge e 

det{D + A)g = ae det{D + A)G-e + det{D + A)G/e- (5.18) 



Proof We pick up a semi-regular edge e entering vi and exiting V2- Thus it exists 

-ei ' ' ej 

•-•-^Xe^^ = l+aeUJeXe. (5.18) 



some i and j with e^l = +1, e^^ = —1. We expand 



leading to two contributions, which we denote respectively by det Qce,! and det QG,e,2- 
For the first term, since one must saturate the Xe and uj^. integrations, one has to 
keep the XeiXvi — Xv2 + Yle feeXe) term and the similar uj term. Note that the sum 
is done on all the edges e hooking to any of the vertices Vi and f 2 and with whom 
the edge e has no trivial Moyal oscillation factor. One has 



detQce,! = y n dxe'dxvduje'dw^ 



e'=^e,v 

iXvi -Xv2+'^ feeXe){^vt -^V2 + ^ fee^e) 

e e 

Q-T,e'j^e'^'eXe"^e' Q-lT,e'^e,vXe'<^c'vXv+X^'^ _ (5-17) 
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As in the commutative case, we now perform the trivial triangular change of variables 
with unit Jacobian: 



Xvi Xvi X' 



V2 



^feeXe, Xv = Xv foT V ^ Vi, (5.17) 



and the same change for the uj variables. What happens now is analogous to the 
commutative case, with the difference that the last term in the definition of 
will lead to the reconstruction of the Moyal oscillation factors of the edges hooking 
to Vi with the edges hooking to V2- This completes the ribbon contraction, thus 
det Qce,! = detQc/e- The second term det(5G,e,2 with the aeOJeXe factor is even 
easier. We must simply put to all terms involving the e label, hence trivially 

det (5G,e,2 = Oie dct Qc-e- ^ 

We need now to compute U* on terminal forms after contracting/deleting all 
semi-regular edges, that is compute on a rosette graph IZ. This is done by using 
the double contraction introduced in the previous section. 

Consider a nice crossing of IZ between two edges ei and 62 with parameters ai 
and It leads to a contribution 

U^ = {aia2 + \e^)Un/e,e, (5.18) 

where we recall that the contracted rosette IZjexe^ is obtained by deleting ei and 
62 from IZ and interchanging the half-edges encompassed by e\ with the ones en- 
compassed by 62, see Figure |9j The procedure can be iterated on 7^/eie2 until after 
gijVj double contractions a planar rosette with 2E{TZ) — 2g{TZ) is reached, for which 
F = and for which the terminal form is Ue as in the commutative case. 

Remark that the main difference with the commutative case is the inclusion of 
the 9^ term in the terminal form evaluation (5.18). This type of genus-term has no 
analog in the commutative case. 



Example 5.1 Consider the graph of Figure 12. Its first Symanzik polynomial is |^ 

aia2 + ttitts + «2a3 + ^^'^^ (5.19) 

Choosing as a regular edge leads to a contracted graph where the pair of edges a\ 
and a2 realizes a nice crossing. We thus have a contribution to the first polynomial 

aia2 + ^9^. (5.20) 
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Figure 12: An example of a non-planar graph, g = 1. 



The deleted part then follows as in the commutative case leading to a contribution 

a^ai + asa2. (5.21) 



Putting together (5.20) and (5.21) leads to the expected result (5.19). 



Let us now give the following definition: 

Definition 5.1 A -k-tree of a connected graph G is a subset of edges with one bound- 
ary. 

This definition allows to write a -k-iiee in some graph of genus g as an ordinary tree 
plus at most g pairs of "genus edges" (where by "genus edges" we understand pairs of 
edges which make a recursive succession of nice crossings under double contractions 
on the rosette obtained after contracting the edges of the tree in the graph). 



Example 5.2 For the graph of Figure 12, the -k-trees are the ordinary trees {1}, {2}, 
{3} and the tree plus one pair of genus edges, namely {1,2,3} which is the whole 
graph. 

In [19], the following general expression for the first polynomial U of the "naive" 
noncommutative model was given 



where we have denoted by 



, as) 



T* it— tree e^T* 



2^ 



(5.22) 



(5.23) 



the number of loops of G. Note that the factor 2 above is the one which matches 
our conventions. 



Let us now give a proof of the formula (5.22). Consider the following lemma: 
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Lemma 5.1 (Lemma III. 2 of j20] ) Let D = be diagonal and 

A = (aij)jjg{i,...,D} be such that an = 0. Then 

det{D + A)= det{Bj^)Y[ai (5.24) 

KC{1,...,N} i(^K 

where is the matrix obtained from A by deleting the lines and columns with indices 
in K. 



The particular form (5.14) of the matrix Q allows thus to use this Lemma to 
calculate its determinant {i.e. the polynomial U). Factoring out | on the first E 
lines and then | on the last V — 1 columns, one has 

^ ^ KC{1,...,E} e(^K 

where we have used that 

b- E = -{V -1). 

Note that the set K on which one sums up corresponds to a set of edges of the graph; 
this comes from the fact that the last V — 1 entries on the diagonal of the matrix A 
are equal to 0. In [20] (see Lemma III. 4) it is proven, using a non trivial triangular 
change of Grassmanian variables that a determinant of type is not vanishing 
if and only if it corresponds to a graph with only one face. This means that the 
complement of the subset of edges K must be a tIt— tree, K = T* . Furthermore, one 
has 

J]ae= n (^-26) 
5.4.2 The translation-invariant model 



The relation of the parametric representation of the model (5.6) to the Bollobas- 
Riordan polynomials follows the one of the "naive" model (5.4) presented above. 
This is an immediate consequence of the intimate relationship between the parametric 
representation of these two noncommutative models, a relationship explained in the 
previous subsection. 
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5.5 The second polynomial for NCQFT 



In this section we prove the form of the second polynomial for the model ( 5.4[ ) (both 



its real and imaginary part, as we will see in the sequel). We then relate this second 
polynomial to the Bollobas-Riordan polynomial. 



From (|5.13|) it follows directly that 

-PRQ-'R'P' (5.27) 



where we have left aside the matrix M coupling the external momenta to themselves. 
Note that the matrix R couples the external momenta (and the hyperposition asso- 
ciated to the root vertex) to the internal momenta and the remaining V — 1 hyper- 
positions. This coupling is done in an analogous way to the coupling of the internal 
momenta with the respective variables. 

We can thus state that the V polynomial is given, as in the commutative case, 
by the inverse Q^^ of the matrix Q giving the U polynomial. The particular form 



(5.14) of the matrix Q leads to 



g-i = l((D + A)-i + (D-A)-i)®i, 

+ ^ {{D + A)-^ - {D - A)-^) ^ a (5.27) 

Thus, the polynomial V has both a real and an imaginary part . 

In the case of the commutative theories, the imaginary part above disappears. 
This is a consequence of the fact that the matrix F, coupling through the Moyal 
oscillations the internal momenta, vanishes for 6 = 0. 

Let the following definitions. 

Definition 5.2 A two -k-tree is a subset of edges with two boundaries. 

Furthermore, let K a subset of lines of the antisymmetric matrix A. Let Vi^Af.-) 
be the Pfaffian of the antisymmetric matrix obtained from A by deleting the edges 
in the set K U {r} for r ^ K. We also define eK,T to be the signature of the 
permutation obtained from (1, . . . , i?) by extracting the positions belonging to K U 
{r} and replacing them at the end in the order: 

> 1, . . . , Z]^, . . . , Zp, . . . , iri • • • 1 if, ip., . . . , i^. 

We now prove a general form for both the real and the imaginary part of the 
polynomial V*, noted and y*. 
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Theorem 5.2 The real part of the second Symanzik polynomial writes 



^ Tj* tree e^Tj* 

where pr* is the sum of the momenta entering one of the two faces of the 2 -k—tree 

■'■2 ■ 

Note that by momentum conservation, the choice of the face in the above theorem is 
irrelevant. Furthermore, let us emphasize on the fact that, being on the submanifold 



Pg = 0, an equivalent writing of (5.28) is 



'^'^ = -2(2)''^' E n2yP.-^^'' (5-29) 

^ ^ r* 2-* tree e^T* 

where p^ (and resp. p^r) is the total momenta entering one of the two faces of the 
2-ic tree. 

Proof. We base our proof on the following lemma: 

Lemma 5.2 (Lemma IV. 1 of |20| ) 

The real part of the polynomial V* writes 



E n f E E Re.rSKr'PKAK,)] (5-30) 



where di are the elements on the diagonal of the matrix Q. Furthermore, when 
\K\ G {E — the matrix with deleted lines is taken to he the empty matrix, with 

unit Pfaffian. 

Note that, as before, since the matrix Q has vanishing entries on the diagonal for the 
last V — 1 entries the subsets K are nothing but subsets of edges. The empty matrix 
obtained from deleting all the first E edges in the graph corresponds to the graph 
with no internal edges but only disconnected vertices. Each of these disconnected 
components has one boundary; hence the Pfaffian is non-vanishing. 



Note that the Pfaffian in (5.30) disappears iff the corresponding graph has 1 



boundary (see above). This means that K U {r} is the complement of a i^r— tree T*: 



KVJ{t] = T\ (5.31) 
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Hence the subset K is the complement of a ^ir-tree plus an edge (just like in the 
commutative case). Adding an extra edge to a ^^r-tree represents an increase of the 
number of boundaries by one unit. Hence, the subset of edges K above is the 
complement of some 2 t^t— tree 7^* 

K = r*. (5.32) 

As before, one has 

The diagonal terms in the matrix Q are again the parameters ae- Factoring out | 
factors on the lines of the matrices corresponding to the edges of the graph and then 
I for the lines of the matrices corresponding to the vertices. The extra factor 9/2 
corresponds to the extra edge r. 

Let us now investigate the square root of the momenta combination entering 



(5.30). Note that the matrix element Re^r is not vanishing only for external momen- 
tum which has a Moyal oscillation with the internal momenta associated to the 
edge r. It is this edge r which actually creates the extra boundary. Thus the sum on 



the external momenta in (5.30) is nothing but the sum of the momenta entering one 



of the two boundaries. By a direct verification, one can explicitly check the signs of 



the respective momenta in (5.30), which concludes the proof. □ 



Example 5.3 For the graph of Figure 12, the second polynomial is 



V*{a,p) = aia2a^p^ + ^(ai + 0:2 + a^)6'^p'^. (5.33) 



Let us now investigate the form of the imaginary part y*. One has the following 
theorem: 

Theorem 5.3 The imaginary part of the second Symanzik polynomial writes 

y*{a,p)=(^-X n2yV^(p), (5-34) 

^ ^ r* * tree e^T* 

where ip{p) is the phase obtained by following the momenta entering the face of the 
•k—tree T* as if it was a Moyal vertex. 
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Proof. The proof follows closely the one of Theorem 5.2 Nevertheless, the equivalent 
of (5.30) is now (see again [20] ) 



^62 



K e^K \ei,e2 \ T,r' / 

(5.34) 

where di are the elements on the diagonal of the matrix Q. Since we look this 
time for sets such that Pf(i?^) is non- vanishing, this implies as above that K is the 
complement of some ^ir-tree T*. Furthermore one needs to consider the two extra 
edges r and r'. It is possible from the initial ^ir— tree above to erase these two more 
edges such that the Pfaffian PfBj^.~, is non- vanishing. Indeed, if the ^ir-tree is a tree, 
by erasing two more edges of it we obtain a graph with 3 disconnected components, 
each of it with a single boundary; the corresponding Pfaffian will be non-vanishing. 
Summing up on all these possibilities leads to the Moyal oscillations of the external 
momenta (the one which disappears when truncating the graph). If the tree is 
formed by a tree and some pair of genus edges we can always delete further the pair of 
genus edges and remain with the regular tree. Obviously the corresponding Pfaffian 
is again non- vanishing (since it corresponds to a graph with only one boundary). □ 
Note that the form of the real part and of the imaginary one of the polynomial 
V* are qualitatively different. Indeed, the real part contains some square of a sum 
of incoming external momenta, while the imaginary one contains a phase involving 
the external momenta. 

Let us end this section by stating that the second noncommutative Symanzik 
polynomial also obeys the deletion-contraction rule. The proof is exactly like in the 



commutative straightforward rereading of Theorem |5.1 



5.6 Relation to multivariate Bollobas-Riordan polynomials 

In the previous subsections, we have identified the first Symanzik polynomial of 
a connected graph in a scalar NCQFT as the first order in w of the multivariate 
Bollobas-Riordan polynomial, 

[/^(a,^) = (^/2)^-^+i(j]a,) X 11111^-1^0(2^,1,^). (5.34) 

eS-B 

Recall that the multivariate Bollobas-Riordan polynomial (see [UH]) is a generaliza- 
tion of the multivariate Tutte polynomial to orientable ribbon graphs defined by the 
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expansion, 



(5.34) 



A<ZE eeA 



with q{A) the number of connected components and b(A) the number of boundaries 
of the spanning graph {V, A). 

In order to deal with the second Symanzik polynomial in the noncommutative 
case, we now introduce an extension of ZG{P,q,w) for ribbon graphs with flags at 
q — 1. In the case of ribbon graphs, the flags are attached to the vertices and the 
cyclic order of flags and half-edges at each vertex matters. For each cyclically oriented 
subset / of the set of labels of the flags, we introduce an independent variable wj. 
Cyclically ordered subsets / arc dcflncd as sequences of different labels up to a cyclic 
permutation. Then, each boundary of a graph with the orientation induced by the 
graph, deflnes a cyclically ordered subset of the set of labels of the flags, by listing 
the flags in the order they appear on the boundary. Accordingly, a variable wj is 
attached to each boundary. 

Definition 5.3 For an orientable ribbon graph G with flags EG{(3e, Wj) is defined by 
the expansion 



where In are the cyclically ordered sets of flags attached to each of the connected 
component of the boundary of the spanning graph (V, ^4) . 

We recover ZG(/3e,l,w) by setting wj = w and ae = 1, but the information 
pertaining to q is lost except for planar graphs. Indeed, in this case the genus of 
any subgraph is still so that \V\ — \A\ + b{A) — 2k{A) and thus ZG{l3e,q,w) — 



The polynomial SG(Q;e, /9e, wi) obeys the contraction/deletion rules for any semi- 
regular edges (i.e. all types of edges except self-loops). The structure of the flags of 
G — e is left unchanged, but less variables wi enter the polynomial since the number 
of boundaries decreases. For G/e, the flags attached to the vertex resulting from the 
contraction are merged respecting the cyclic order of flags and half-edges attached 
to the boundary of the subgraph made of the contracted edge only. 

Proposition 5.1 The polynomial Eaioie, PejWi) obeys the contraction/deletion rule 
for a semi-regular edge, 

EG{ae,Pe,Wi) = aeEo-eicte'^e, Pe'^e,Wl) + Pe'^G/e{ae'^e, Pe'^e, Wl) ■ (5.34) 




(5.34) 



?l^l/2ZG(g-Me,g5«;). 
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This follows from gathering in the expansion of ^(^(ae, /5e, w/) the terms that 
contain e and those that do not. The contraction/deletion rule may be extended to 
any edge provided we introduce vertices that are surfaces with boundaries as in [68] . 

The second interesting property of Eciae, Pe,wi) lies in its invariance under du- 
ality. Recall that for a connected ribbon graph G with flags, its dual G* is defined 
by taking as vertices the boundaries of G, with flags and half-edges attached in the 
cyclic order following the orientation of the boundary induced by that of G. 

Proposition 5.2 For a connected graph G with dual G* , 



Proof: First, recall that there is a natural bijection between the edges of G and 
those of G* . Thus, to a subset A of edges of G we associate a subset A* of edges 
of G* which is the image under the previous bijection of the complementary E — A. 
Then, the term corresponding to A in ^^(ae, /3e, "W^/) equals that corresponding to 
A* in SG*(/5e, cxei wj). The only non trivial part in the last statement is the equality 
of the boundary terms in G and G*, which is best understood by embedding G in a 
surface S. Then, the spanning graph {V*,A*), viewed as discs joined by ribbons, is 
homeomorphic to S — {V,A), with the orientation reversed. Accordingly, they have 
the same boundary. 



This relation may also be extended to non connected graphs at the price of 
introducing again vertices that are surfaces with holes. For example, the dual of 
a disjoint union of n vertices is the vertex made of a sphere with n holes. For a 
regular edge, the duality exchanges contraction (resp. deletion) in G with deletion 
(resp. contraction) in G*. In the case of the deletion of a bridge in G, we have 
to contract a self-loop in G*, thus leading to vertices that are surfaces with holes. 
Note that this implies a duality for the multivariate BoUobas-Riordan polynomial 
only at the special point g = 1, in agreement with the fact that the duality for the 
Bollobas-Riordan polynomial only holds when its arguments lies on a hypersurface 



Finally, let us come to the relation with the second Symanzik polynomial in 
NCQFT. For a given connected graph with momenta pi such that YliPi — attached 
to the flags, we decompose the latter polynomial into real and imaginary part. 



(5.34) 



□ 



m 



'a{ae,9,pi) = X^{ae,9,pi) +iy^{ae,9,pi). 



(5.34) 
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Consider real variables Wi and define wi = Wi for any cyclically oriented subset 
of flags, Then, expand (6'/2)''^'~'^' 'EG{2ae/0, 9wi/2) to the first two orders at Wi = 0, 

(e/2y^\~\''\EG{2aJ9,9wi/2) = Af^J^"^,) B^JW,w, + o{w^y (5.34) 

The first order term reproduces the first Symanzik polynomial 

f/^(ae,^) = A (5.34) 

whereas the second order terms yields the real part of the second Symanzik polyno- 
mial, 

A'*(ae, 0,pi) = -\Y^ AijPi ■ pj. (5.34) 
To obtain the imaginary part, consider the variables 

i<j 

if / contain all the flags and wi = otherwise. The previous definition involves 
a choice of a total order on / compatible with its cyclic structure, but momentum 
conservation ^jPj = implies that does not depend on this choice. Then, 

yUc^e,0,p,) = i9/2y^\-\''\EG{2aJ9,wi). (5.34) 

As a consequence of their expressions in terms of Edae, Pe, wj), the noncommutative 
Symanzik polynomials obey contraction/deletion rules for regular edges and duality 
relations. For example, the duality for the first Symanzik polynomial reads 

{9/2r\ U^aiae, 9) = {9/2r\ ( H ^) (9^1 a., 9) . (5.34) 

eeE 

Beware that G* is the dual graph whereas the star on polynomials such as f/* and V'^ 
refer to the Moyal product. Analogous relations, though slightly more cumbersome, 
can be written for the second Symanzik polynomial. 

Still an other way to categorify and regularize in the infrared is to introduce har- 
monic potentials on the edges rather than the vertices, leading to propagators based 
on the Mehler rather than the heat kernel. This is the so-called vulcanization. An 
extensive study of the corresponding commutative and noncommutative polynomials 
is under way as a companion paper |7Uj . 
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